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Abstract 

Q\ The isomonodromic tau function defined by Jimbo-Miwa-Ueno vanishes on the Malgrange's divisor of generalized 
monodromy data for which a vector bundle is nontrivial, or, which is the same, a certain Riemann-Hilbert problem 
has no solution. In their original work, Jimbo, Miwa, Ueno provided an algebraic construction of its derivatives with 
^ respect to isomonodromic times. However the dependence on the (generalized) monodromy data (i.e. monodromy 
representation and Stokes' parameters) was not derived. We fill the gap by providing a (simpler and more general) 
description in which all the parameters of the problem (monodromy-changing and monodromy-preserving) are dealt 
with at the same level. We thus provide variational formulas for the isomonodromic tau function with respect to 
the (generalized) monodromy data. The construction applies more generally: given any (sufficiently well-behaved) 
family of Riemann-Hilbert problems (RHP) where the jump matrices depend arbitrarily on deformation parameters, 
we can construct a one-form (not necessarily closed) on the deformation space (Malgrange's differential), defined 
off Malgrange's divisor. We then introduce the notion of discrete Schlesinger transformation: it means that we 
allow the solution of the RHP to have poles (or zeros) at prescribed point(s). Even if fl is not closed, its difference 
evaluated along the original solution and the transformed one, is shown to be the logarithmic differential (on the 
deformation space) of a function. As a function of the position of the points of the Schlesinger transformation, 
yields a natural generalization of Sato formula for the Baker-Akhiezer vector even in the absence of a tau function, 

^*>0 and it realizes the solution of the RHP as such BA vector. 

Some exemplifications in the setting of the Painlcve II equation and finite Toplitz/Hankel determinants are 

provided. 
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1 Introduction 

In the eighties Jimbo, Miwa and Ueno [9, 7, 8] derived a set of algebraic (in fact rational) nonlinear equations 
describing deformations of a rational connection on P 1 which preserve the generalized monodromy data. They 
associated to this deformation a closed differential u> JMU on the space of deformation parameters, namely, on the 
space of " isomonodromic times" which we denote collectively by t. In the simplest case of Fuchsian singularities 

K A 

*'(*) = , *(oo) = 1 , A(z) = £;— £)A,- = ( L1 ) 

3=1 Z ° j 

the isomonodromic deformation equations were studied by Schlesinger [15] 

6A k = -Y l [Ak,A j ] S{ak ~ aj) , 5:=J2^jd aj (1.2) 

and the Jimbo-Miwa-Ueno differential reads 

«W = I E Tr(^)fc^ (1.3) 

*.../< a j a k 

3,kj^=k J 

It can be checked directly that if the matrices Ak(a) depend on the position of the poles as mandated by (1.2) then 
uj jmu above is a closed differential. 

This differential was generalized to arbitrary (generic) rational connection in [9], to which we refer the reader 
for details. In the above situation for Schlesinger deformations, the locations of the poles constitute the "isomon- 
odromic" parameters or times and we denote them by t. This was an important achievement because of the sweeping 
applications of isomonodromic deformation to integrable systems (solitons solutions to KP, solutions to Toda, etc.) 
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Painleve equations and, later, random matrix models. The (exponential) integral of this closed differential is the 
"isomonodromic tau function" 

T JMU $m)=eS"JMv . (1.4) 

and the Painleve property translates to the fact that t jmu is a holomorphic function of the isomonodromic times 
that has only zeroes (away from an explicit set of times where it has a branching behavior: in the case of (1.2) 
this is the set of "diagonals" aj = a^, j ' ^ k). In (1.4) we have indicated that the tau function depends necessarily 
on the (generalized) monodromy data, denoted generically by m; this dependence is parametric and the present 
paper addresses precisely the 

Question 1 (Naive) What is the dependence of t jmu on the monodromy data? 

The question is conceptually simple but slightly ill-posed; since what JMU really defined was only a differential 
(in symbolic notation) 

u.,mu = fj (*? fn) d*j , d tj f k =dtjj (1.5) 
j 

the dependence of t on m is defined only up to multiplication by an arbitrary function of the fh's only. So a better 
question would be 

Question 2 (Refined) What is the essential dependence of T JMV on the monodromy data?. Can we define an 
extended closed differential £1 on the total phase space of the problem that coincides with io JMU on the isomon- 
odromic submanifold, namely 

Uext = m)dtj + y^ j G v {t,fh)dm v (1.6) 

d m „fk = 9t k Gv , dm^G^ = dm^Gfj,. (1.7) 



The question still admits many solutions as stated, since if we find one such extension we are still free to add any 
closed differential of the m's alone. However we may and should understand the problem in a relative setting, where 
the answer is taken modulo closed forms of the m's which must be holomorphic on the whole space of monodromy 
data. The reason for requiring holomorphicity is actually important because of the interpretation of the singularity 
locus of u>, as we presently explain. 

The meaning of the singularity locus of lu jmu . Malgrange (for Fucshian systems) [11] and later Palmer 
(for irregular singularities) [14] showed what the meaning of the zero-locus of t jmu is: when T JMU (t; rn) — then 
a vector bundle on P 1 is nontrivial or -which is equivalent- a Riemann-Hilbert problem is not solvable. This 
is the equivalent of saying that lo JMU has only simple poles (with "residue" one) away from the non-movable 
singularity locus (Painleve property). The divisor where the aforementioned bundle is nontrivial is generally 
termed Malgrange divisor (or simply Malgrange divisor). 

It is clear then that whatever extension (1.6) we are looking for, it ought to preserve the singularity locus, that 
is, the G v (t, in) may have singularities only where some of the fk's has, and of the same type. 

So we arrive to the final formulation of a "sensible" problem, whose solution is the principal aim of the paper 
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Problem 1.1 Formulate a ("natural") extended closed differential oj ext (1-6) such that its tau function (locally 
defined up to nonzero multiplicative constants) 

T ext (t,fh)=eJ u °** (1.8) 

vanishes precisely and only on the Malgrange divisor in the extended phase space of isomonodromic times t and 
monodromy data fh. 

The differential we propose is in fact very natural (see Def. 2.2 and Thm. 5.1): it is the pull-back of the (partial 
integral of the) generator of the third cohomology of the "loop group" to the submanifold corresponding to our 
total phase space. In particular it expresses the same cohomology class c[0] [11]. 

We will briefly indicate some interesting problems which can be addressed and that require the knowledge of 
the derivative of t with respect to non-isomonodromic times (hence the knowledge of the G^'s). 

Organization of the paper. The heart of the paper is in fact Section 2 where we introduce (recall) the definition 
of the Malgrange differential lu m (Def. 2.1) associated to any (sufficiently well-behaved) Riemann-Hilbert problem. 
Such differential is ill-defined only when the RHP is not solvable: however its exterior differential is (admits a) 
smooth (extension) over the whole "loop group" , in particular also at the points where the RHP is not solvable 
[11]. The curvature of lo m is not zero, but it is so explicit that it is immediate to identify "simple" families of 
Riemann-Hilbert problems for which ui M is closed. By adding to oj m an explicit smooth differential (in particular 
this does not change its singularity locus) we obtain a new differential (Def. 2.2) whose curvature differs slightly 
from that of lo m but within the same cohomology class. 

In Section 2.3 we investigate the changes in f2 (or uj m ) under modification of the growth conditions for the 
solution of the RHP (discrete "Schlesinger" transformations); this allows to interpret the matrix solution of any 
Ricmann Hilbcrt problem as a Baker-Akhiczcr function via a Sato-like formula, even if a notion of tau-function 
is not available (see Thm. 2.1 and Sect. 2.3.1, Sect. 2.3.2). 

In the second part we specialize the setting from arbitrary Riemann-Hilbert problems to those that correspond 
to rational ODEs in the complex plane: Section 3 is a quick reminder to the reader of the classical description 
of the (generalized) monodromy map, i.e. how to associate to a rational ODE in the complex plane the set of 
"Birkhoff data" of irregular type, connection matrices, monodromy representation and Stokes' multipliers. All the 
material is quite standard [17]. 

These Birkhoff data can be used viceversa (Sect. 4) to encode an ODE in a Riemann-Hilbert problem (which 
may or may not have a solution, although generically it does [16, 17]). 

In Section 5 we show that f2 is a closed differential in all the deformation parameters, which include 
• Monodromy representation; • Connection matrices; • Stokes' matrices. 

In Section 5.1 we show that the restriction of f2 to the submanifold of isomonodromic times coincides with the 
differential oj jmu defined in [9]. Therefore £1 realizes the solution to our Problem 1.1. 

We conclude the paper with Section 6, where applications are provided to Painleve II, (shifted) Toeplitz and 
Hankel determinants. 

Acknowledgements. The author thanks Jacques Hurtubise and John Harnad for helpful discussion. The work 
was partially supported by the Natural Sciences and Engineering Research Council of Canada (NSERC) . 
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2 Riemann— Hilbert problems 



A Riemann-Hilbert problem starts with the following data. We will assume latitude in the smoothness class as 
this is not our primary focus. 

The Riemann Hilbert data 

1. a finite collection of smooth oriented arcs j„, j — 1 . . . K, possibly meeting at a finite number of points but 
always in non-tangential way (Fig. 3 is a good example). We denote collectively these arcs by the symbol £7 

2. a collection of r x r matrices Af„(z), each of which holomorphic in a neighborhood of each interior point of the 
corresponding arc 7„. We make the assumption that these matrices have unit determinant 3 det M v (z) = 1. 
We will denote collectively by M(z) the matrix defined on £7 that coincides with M u (z) on j U) 

M : £7 -> SL r (C) 

z ~ X)M„(z)x 7 „(z) (2.2) 

where, for a set S, \S denotes its indicator function. 
The Riemann Hilbert problem then consists in finding a holomorphic matrix on the complement of the contours 

T(z) : C\£ 7 -> SL r (C) (2.3) 
such that it admits (non-tangential) boundary values satisfying 

T + (z) = T_(z)M(z) ze£ 7 . (2.4) 
The above data are insufficient to characterize uniquely the matrix T (if it exists) and need to be supplemented by 

• growth behavior near the endpoints/intersections of the contours 7„ and at co; 

• an overall normalization. 
Therefore we will pose the following 

Problem 2.1 (RHP) Find a holomorphic matrix T : C \ £7 — > GL n (C) such that 

• T+(x) = T-(x)M(x) x e £7; 

• T(z) is uniformly bounded in C; 

• T(zo) = 1 

Typically we will choose zq = 00. 

3 This is not really necessary but simplifies some matters. In some situations (typically where the determinant a rational function) 
it is necessary (but simple) to relax the condition and allow any matrices in GL r (C). We face the problem on a need-to basis. 
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In general there are obstructions to the solvability of 
Problem 2.1: for example let c be a point in £7 where 
several arcs 71, ... ,7^ meet (at nonzero tangential rel- 
ative angles) and suppose they are oriented towards c. 
Let Mj (z), j — 1 , . . . , £ are the corresponding (locally 
analytic) jump matrices, then 

M 1 {c)M 2 {c)...M e {c) ? 1 . (2.5) 

is an obstruction. This does not mean that it is im- 
possible to find r satisfying the jumps, but it cannot 
be bounded at c in general and thus a relaxed growth 
constraint must be allowed at c. 

However there are situations which are of interest 
for our applications where a growth behavior can be 
"traded in" for an extra contour. 

Example 2.1 A typical example is depicted in Fig. 1; here on the contour that terminates at z = a there is a constant 
jump of the form M — C~ 1 e 2tnL C , where L is upper triangular (in Jordan form, for example). It is clear that the diagonal 
entries of L (the exponents of formal monodromyj are defined only up to addition of integers. This arbitrariness in 
fact corresponds to the choice of a growth behavior at the endpoint for the solution of the RHP. Once this choice has been 
made, we can always recast the problem into an equivalent one where the solution V is bounded. This is achieved by adding 
a small circle at a and re-defining F — FC~ 1 (z — a)~ L inside: in the new formulation for the problem, F will be bounded in 
a neighborhood of z = a and the two problems are equivalent, in the sense that one has solution if and only if the other does, 
and the relation between the two solutions is also very simple. 

(Loop) group structure. Although we will not make any explicit use of the following fact, we mention that the 
matrices M : £7 — > SL r (C) that satisfy all the condition above form a (infinite dimensional) Lie group, akin to 
the usual loop group, were it not for the fact that the contour is not a (collection of) circle(s). It is convenient to 
introduce a symbol for this manifold (group) 

G := j M : U 7l/ h-> SL r (C), (z) locally analytic j (2.6) 

Since we will only consider finite dimensional (analytic) submanifolds of Q we will not dwell on the infinite- 
dimensional differential-geometric issues. Most notably the tangent bundle TQ will be used in a rather naive form 
where possible issues of topological nature (arising from the infinite dimensionality) are disposed of. In fact only 
the tangent space to the finite dimensional submanifolds of interest will appear, and hence the point is not relevant. 

2.1 Deformations and Malgrange's form 

Suppose now that M(z) = M(z; s) depends holomorphically on additional parameters s; the reader should think of 
this as an explicit dependence, dictated by the problem under consideration. The parameters s could be thought 
of as coordinates on a manifold of deformations. When -eventually- we specialize the setting these s's will be the 
isomonodromic times together with the monodromy data. On this manifold we define the one-form (differential) 
already used by Malgrange [11]. 



1. 

T~Y(z)(z-a) L C 




Figure 1: The process of "cloaking" of a point of growth. 
The two RHPs are equivalent in the sense that one admits 
a solution if and only if the other does. 
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Definition 2.1 Let d denote the derivative w.r.t. one of the parameters s. Then we define Malgrange's form lu m 

u M (d)=LJ M (d;[T}):=-^Tr(rz\x)rL(x)E d (x)^ (2.7) 

Eg(z) := dM(z)M~ 1 {z) . 

Remark 2.1 The minus sign is -of course- conventional, but it is important because with the minus the tau- 
f unction (when it exists) has a zero and not a pole on the Malgrange divisor. 

Remark 2.2 The definition would hold identically also for any RHP formulated on a Riemann surface. 

Remark 2.3 The reader may frown upon the boundary value of the derivative ofT: however, since our assumptions 
mandate that M v (z) is analytic in a neighborhood of the contour 7„, it is simple to show that the derivative V also 
admits (bounded) boundary values. 

Remark 2.4 In the second notation for lo m we have indicated the dependence on T, which depends on the choice 
of jump matrices and specifications of the growth behaviors. When no ambiguity can arise, we will understand 
such dependence. 

Curvature of uj m . The first issue is the computation of the exterior differential: the computation can be found 
in [il] but in rather abstract terms and we prefer to give a direct derivation here. We have first the 

Lemma 2.1 Let d denote any vector field in the parameters of the jump matrices. Then 

dr(z)=[ r-WWW ^^) (2 . 8) 



\j- X — Z 2«7T 



and, consequently 



8 (^EM) = f r-H^(x) Ed{x) Tz\xn w) ax_ 



z — w ) ,/s 7 z — x x — w 2in 



Proof First of all we have 

dT + = dV-M + T-dM . (2.10) 

This is a non-homogeneous Riemann-Hilbert problem: for convenience we take the normalization point at infinity 
r(oo) = 1 so that dT(z) = 0(z~ 1 ). It is promptly seen that the proposed expression fulfills (2.10) and this last 
condition. The uniqueness follows from the fact that the homogeneous part has only the trivial solution that tends 
to zero at oo. The formula (2.9) follows by direct application of (2.8) and Leibnitz rule. Q.E.D. 

Proposition 2.1 Denoting with d,d the derivatives w.r.t to two of the parameters s, the exterior differential of 
uj m is 

r,(d, d) := dco M 0) - dco M (d) =\J^ (%(x)-^~ 9 (x) - ^-~s(x)~ d (x)^ ^ (2.11) 
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Proof Let d, d be two commuting vector fields so that 

8u(d,d) = du{d)-dw(d). 

We have 

da; 



du M {d) = -d / Tr (rz\x)T'_(x)E d (x)) - 



E 7 



2m 



Tr(d(TZ\x)T'_(x))Eg(x))— + Tr(TZ 1 (x)T'_(x)d(E d (x))) 



Tr 



E 7 



dy 



da: 
2i?r 



2in 

y)-d(y)£(y, xJE^x) 

E 7 (x~y) 2 j x=x± ^n JSl 

where we have introduced the convenience notation 

}C(x,y) = TZ 1 (x)T_(y) . 

Continuing our computation we find 



2m 



(2.12) 

(2.13) 
(2.14) 



Tr 



{Yz\x)Y'_{x) (d3 a {x)+[3 B ,E s ])ty 



2m 



(2.16) 



du u (d) - du M {d) = 



Tr 



E 7 JE 7 



n ( IC(x,y)Eg(y)K.(y,x)E e (x) ^ 



Tr 



£7 </£7 



(x - y) 
y)E d (y)IC(y, x)E§(x) 



E 7 



Tr I I 1-1 ! 1 ' 



(x - y) 2 
Ax 
2in 



dy 



dx 

2in + 
dx 

. 2i7T^ 



(2.17) 



In the first two terms the order of integration is important since the kernel is singular due to the denominator 
(x — y) 2 . Note that -by our assumptions- the jump matrices have analytic continuation in a neighborhood of the 
contours. A standard computation involving residues yields the expression 



Tr 



^0, y)-s(y)fc(y, x)E d {x) 



E 7 JE7 



(x - y)' 



dx 



dy 



Tr fC(x, y)E d (y)fC(y, x)E s (x) 



) ^2m 
dx 1 
,2iTT ~ 2 



Tr 

E 7 ^S 7 V 

Tr( E a (x)Eg(x) 



K>{x, y)E a {y)K{y, x)E s (x) 



{x- y y 



dy 



da; 

27tt ~ 
(2.18) 



Tr [TZ 1 (x)T'_(x)E d (x)E B (x) - Ee^TZ 1 (x)T>_(x)E B (x) + E' a (x)E B (x)] ^ - ^Tr (s a (a;)%(a;)) 



dx 
2m 

Tr(TZ l {x)T'_{x) [E d (x),E d (x)] ) ^ + \ jf Tr (e' 9 (x)E b (x) - E a (x)E' B (x)) 



E 7 



Tr (rz\x)r'_(x) [E d (x),Eg(x)}) ^ + / Tr (S' a (a;)%(a;)) g - ±Tr[ E„(.,-.E ,l , I 



a;e5E 7 



da: 

2i7T 



3tE<5£7 

(2.19) 
(2.20) 



The notation for the last term in (2.18) is the evaluation of the "boundary term" at the endpoints of all the 7„ 
according to the rules of integration by parts (i.e. with the appropriate signs depending on the orientations of the 
contours) . 

Thus, resuming (2.17), we have 



r](d,d) := du M {d) - duj M {d) 



s 7 



Tr(S 5 (x)— E 9 {x)- 



Eg(x)E a (x) 



dx 

2i7T 



(2.21) 



Q.E.D. 



s 



The Theta divisor. The reader should observe and compare the expressions (2.11) and (2.7): the crucial point 
is that ojm (2.7) is not defined whenever the RHP does not admit a solution. On the contrary, as a sort of "miracle" 
the expression (2.11) for its curvature is defined and holomorphic for any group-valued matrix M(z), whether or 
not the RHP is solvable. In other words, r\ is a closed two-form which is smooth also on the Malgrange divisor, 
whereas ujm is undefined at those points. Malgrange proved [11] (and Palmer generalized [14]) that w« has a 
"pole" on such divisor: evidently, in the curvature such a pole disappears. 

Tau function It is also apparent that ujm is in general not a closed one-form on the whole (infinite dimensional) 
manifold Q: however, it may become closed when restricted to suitable submanifolds C > Q. Since the curvature 
rj (2.11) is explicitly computable, it is easy to verify for a given explicit submanifold C whether u>m \ c is closed or 
not. When such restriction turns out to be closed we can define (locally) a function on C by 

T C = e JuJ "\r. (2.22) 

By Malgrange's (Palmer's) results, this function (defined up to nonzero multiplicative constant) will vanish at the 
intersection with the O divisor 

t c (s) = o & seCne (2.23) 

Suppose, however that to M | is not closed; we may still seek another differential 9 that "cures" the curvature 

d(o; M +0)| £ = O (2.24) 

and then proceed with the construction of the tau function as before. 

Ideally such differential should be smooth on the whole Q so that it does not change the cohomology class of 
ujm and also does not change the singularity structure; this way the tau function will still vanish only at C n and 
our goal is met. 

Although it is not a pedagogical approach, since we know what C will be for our purposes and we have already 
found 0, we will describe it directly here. 

Changing the curvature of uj m . Consider the one form that, evaluated on a vector d yields 

1 t dr 

0(0) : = — / Tr (M'M dM M ) — (2.25) 

2 jst 2i7r 



Its curvature is the two-form 



60(d, 0) = 00(0) - d9(d) = - [ Tr (dM'M^dMM- 1 - DM' M~ 1 dMM~ 1 ") — (2.26) 

2 Js 7 V / 2m 

which the reader can verify using Leibnitz rule; here the prime denotes (as always) the derivative with respect to z. 
The important (but trivial) additional observation is that is a smooth differential on the whole G, since it does 
not require the solution of a RHP. 

Definition 2.2 The modified Malgrange differential is defined as CI :— lj m + 0. 
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Remark 2.5 A direct computation shows that 

f Tr (r^r'JMM" 1 ) — - / Tr (T^T'+M^dM) — = - f Tr (M'tf^SMM" 1 ) — = 20(0) (2.27) 

so f/iai f2 can be written in the more symmetric form 

i r C ] T 
n(d; [T]) = -- / Tr (TZ 1 T'_dMAl~ 1 + T+ 1 T' + M~ 1 dM) — (2.28) 

Combining the computation of the curvature of u) M (Prop. 2.1) with the curvature of $ above we have the simple 
Proposition 2.2 The curvature of the modified Malgrange form is 

Sn(d,d) = 00(0) - da(d) = - [ Tt( WM- 1 \dMM~ 1 , dMM -1 ] ) — (2.29) 

2 Je 7 V L J / 2«7r 

Remark 2.6 Referring to [11] Thm. 5.5 the curvature 5£l clearly lies in the same cohomology class as 5u> M since 
the two differ by a smooth exact differential. If we define 7 as the total Maurer-Cartan form 7 = M' ' M~ x dx + 
J>2 djMM~ x dsj one has the expression 



5Q = f Tr( 7 A 7 A 7) (2.30) 

127TI ./e 7 



2.2 Submanifolds of Q where f2 is closed 

Looking at the formula of the curvature form of Prop. 2.2 a nontrivial class consists of M(z) that on each arc of 
£7 reduce to one of the following forms: 

• (Piecewise triangular) Matrices of the forms 

M v (z) = P„(l + N^P- 1 , z6> (2.31) 

where N v (z) are upper-triangular analytic matrices and P v is any constant permutation matrix (i.e. an 
element of the Weyl group for SL r ) 

• (Constants) matrices independent of z; 

• (Torals): matrices M„ = D u (z) with D v (z) diagonal matrices or any conjugation thereof by an arbitrarily 
chosen but fixed matrix. 

We see in Section 3 that any (generic) rational ODE can be encoded in a Riemann-Hilbert problem with jumps 
of the form indicated here above. Therefore il on these submanifolds yields a closed differential. 
We will also show that its restriction to the isomonodromic submanifolds coincides with uj ...... 
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2.3 "Schlesinger" transformations 

The aim of this section is to compare the differential cu M on two RHP for r,T defined by the same jumps but 
different growth behaviors, in the form of integer powers for columns at various points. We note immediately 
that two such solution differ (if both existent) by a left multiplication by a rational matrix R(z); indeed -having 
r, r the same jumps by assumption- the matrix 

R{z) :=f(z)r- 1 (z) (2.32) 

is an analytic function in C taken away the points where T has different growth. If this difference in growth is 
polynomial (as we seek now) then R(z) is forced to be rational. 

We now make the following observation: let R(z) be a rational matrix such that the divisor of all poles of 
R, R~ x consists of the point ci, . . . , cr- ^ £7. Define 

f(z) := R(z)T(z) (2.33) 

Quite clearly T solve a different RHP with the same jumps but different growth at the poles of R. The difference 
between the 10 M evaluated along the two different solutions is given by 

uj m (d; [f ]) - u u (d; [T]) = - [ Tr (f I 1 ? 'JMM" 1 ) + f Tr (r^r'JMM" 1 ) = 

J £7 ^ ' J S7 

Tr (TZ 1 R~ 1 R'T_dMM- 1 ) = - [ Tr (R^R'T^dMM^TZ 1 ) = 
E7 J S7 

K 

= / Tr (r^r' (ar+r; 1 - <9r_r_)) = V res Tr (nr 1 marr- 1 ) (2.34) 

JSy j = l Z ~ Cj 

Note that -since fl (Def. 2.2) differ by u) M only in an explicit term d that depends only on the jump matrices, we 
have 

cu M (d; [f]) - u M (d; [T]) = Q(d; [f]) - n(d; [T]) (2.35) 
We now proceed with the definitions in the title of the section and specialize the class of rational matrices R(z) . 
Definition 2.3 Given two distinct points £ ^ r\ and two (possibly formal) series 

(OO \ / OO \ 

1 + , Y v {z) = G n U + ^Y me zA (2.36) 

z x := (z - x), Zoo := - (2.37) 

z 

an elementary Schlesinger transformations at two distinct points £ 7^ 77 is a rational matrix R(z) such that 

R{z)Y i {z) - %{z)zf , R(z)Y v (z) = %{z)z^ (2.38) 

where Y, denote formal series of the same form. In the case of 00 we have Goo = 1 = Goo. If neither £ =/= 00 =/= r\ 
then we impose also R(po) = 1. 

[For the case £ = rj the definition should be modified in an obvious way: please see below]. 

The problem is purely of algebraic nature, and not a very difficult one: the computation is contained in [[7], 
App. A] and the derivation will not be reported (we will give below the relevant results). 

Suppose now we have a RHP for T with jump matrices M as in Prob. 2.1: at points £ ^ £7 the solution T to 
Prob. 2.1 yields a (convergent) power series which can be used as input for the above procedures. 

If £ G £7 we cannot have a (even formal) series since -in general- not even the value of T is well defined at z = £. 
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In the applications to ODEs the following situation occurs; let £ S £7 be 
a point where 71, ... ,7^ meet > 1). Suppose that the jump matrices are 
analytic in small sectors centered at £ containing the direction of approach 
and that -in said sector- (see Fig. 2) 



72 



M 3 {z) = l + 0{{z-tr) 



(2.39) 



Then it is not hard to see 4 that the solution T(z) has the same asymptotic 
expansion in each of the sectors at £ separated by the incoming arcs 



74 



3 — 



(2.40) 



Figure 2: Illustration of the sectors 
of analyticity for Mj where the de- 
cay (2.39) should be valid. 

We will allow to perform Schlesinger transformations involving either points 

£ not in £7 or points where the condition (2.39) is met, so that the solution T defines unambiguously a (formal) 
analytic series centered at the point. 

Definition 2.4 The elementary Schlesinger transformation 1 1 7 j for the solution T of Problem 2.1 is the solution 
(if it exists) of the following new RHP (where z^ = (z — £) if £ 00 and z^ := -) 







t = v 


(* * J) 


?+ = 

f(z) 
f(z) 


= r_M 

= 0(l)z v E i , z~ v 


r+ = 

f(z) 


r_M 

= 0(l)z ? 


-E i+Ej ; z „ £ 



(2.41) 



where the normalization is fixed by requiring that T(z) ~ 1 + 0(2; 1 ) i/£ 7^ 00 7^ 7/ or -if either £ or 77 are infinity- 
that the term indicated by 0(1) above is actually 1 + C(z _1 ). 

Here and below we use the notation Eij for the elementary matrices (with a 1 on the i-th row, j-ih column) and 
Ei := En. It is immediately seen that 

f(z) = R(z)T(z) (2.42) 

with R{z) rational: indeed the ratio T(z)T^ 1 (z) for the transform | ^ 7 j- does not have jumps and may have at 
most a simple pole at z = £. The matrix i?(z) can be computed directly from the (possibly formal) series-expansion 
of r at £,77 (see [7]. We give below a more compact formula that the reader can check autonomously 

4 By moving slightly the jump -which can be done due to the local analyticity of the jump matrix— one sees that the solution Y can 
be "continued" analytically across the jump from the left and from the right. The ratio of these two extension in the common sector 
of analyticity differs from the identity by exponentially small terms, which are transparent to any asymptotic expansion. 
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Proposition 2.3 (Cf. App. A, B in [7]) The left-multiplier matrix R{z) that implements the elementary Schlesinger 
transform 1 1 ? j is given by 



i?(z) = 1 + 

fr x (z) = i 



(r-i(Or(r?))y *-£ 

(r-HOr^))^ z-rj 

z — rj 



det R(z) = 



z~Z 



(ill jgj) 



= 1 + 

ir 1 ^) = i 



l r(a^r-Hg) 
(r-i(Or'(O)ti 

i r(0^r- 1 K) 



(ill (j^j 



i?(z) = 1-E j 



Ri 



Ri 



-^ n [(i-E jj )r'(oo)-i]E ji r- 1 {0 



R- 1 {z) = E n {z-£) + R 



det i?(z) = 



(r-HOr'(O)u 

det .R = 1 



R{z) = E u (z - 7]) + R 
R = {l-E ii T'{ OO )) (l ' ; 



i?- 1 (z) = l-^ + ^ 

i?i = r^-r^)^ [r'(oo)(i - + 1] 

det R(z) = (z — rf) 



where we have denoted by r'(oo) the derivative in the local parameter, that is 

Viz) = 1 + -r'(oo) + . . . . 

z 

(the formulas for | °° °° | can oe found in loc cit.) 



(2.43) 



(2.44) 



(2.45) 



If £ (or r\ or both) belong to £7 where condition (2.39) is in place, the computation leading to (2.34) can be 
still carried out with minor modifications in the process but not in the result. First of all note that the integral 
0J M (d, [r]) is still convergent since Eq = dMM~ Y — 0((x — c)°°) (recall that we assumed M = l + 0((x — c)°°)) and 
hence the additional algebraic growth of r _1 r' along the contours incident at c is still integrable when multiplied 
by Sa. 

Let S7 e denote the contours that lie outside of e disks centered at the points z = £, r) (possibly the same) of 
the Schlesinger transform; we then have 

u M (d; [f]) - u> M (d; [T]) = Km J Tr ((f I 1 ? - I^Pl) Eg) ^ = 
= - lim / Tr (TZ 1 R^R'T-dMM- 1 ) = - lim / Tr {R~ 1 R'T- dMM~ X YZ X ) = 
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= Km / Tr (R^R 1 (dT+T^ 1 - <9r_r_)) = Km I Tr (R^R'dTT- 1 ) (2.46) 

£ ^°Js 7s e ~> J\z-Z\=\z-T)\=e 

In the computation of the limit (2.46) we can replace cffT -1 by a suitable truncation of the asymptotic series at 
z = £,77 (which do not depend on the direction of approach under our assumptions for M), committing an o(l) 
error as e — > 0. Then the limit equals the formal residue 

oj M (d; [f]) - oj M (d; [r]) = res Tr (r^R'OY^ {z)% (z))" 1 ) dz + res Tr [r- 1 R' dY^Y^z))- 1 ^ dz (2.47) 

Here the residue means simply the coefficient of the power —1 of the local parameter (note that has at most 

a double pole and hence the formal residue involves at most the first two terms in the formal series Y). In case 
only one point is involved in the Schlesinger transformation we have simply only one residue at the end. 

Remark 2.7 In the case of a Schlesinger transformation involving two distinct points the determinant of the 
solution det T cannot remain constant, since R(z) has non-constant determinant. This does not pose any significant 
problem as we explain presently. 

The modification can be explained as follows: to account for the different power-law of the columns at the two 
points z — £, z — r\, small counterclockwise circles around those points should be added to £7 imposing additional 
jumps of the form 

M^z) = {z - £)- Ei , \z-£\=e, M 2 {z) = {z-rj) E * , \z-r,\ = e, (2.48) 

conjugating by the same matrices any jump of a contour that passes within said circles. Of course these jumps do 
not have unit determinant and hence the uniqueness of the solution must be argued in a different way from the one 
used in the unimodular case. However the modification in the reasoning is only minor: for a solution T(z) of the 
new RHP we see that detr(z) is analytic and bounded everywhere, except for jumps on the new circles where 

detr+(z) = detr_(z)(z - £) _1 V - CI = e and det 1 + (z) = det T_ (z) (z - if) , \z-i)\ = t (2.49) 

This means that detr(z) admits analytic continuation in the interior of the two disks, with a simple pole at z = n 
and a simple zero at z — £, plus the condition detr(oo) = 1. This forces detT(z) = outside of the disks, 

detr(z) = {z — £) for \z — T)\ < e and viceversa detT(z) = for \z — £| < e. Any solution will have the same 
determinant and hence the uniqueness is established along the same way used previously. 

A direct computation based on (2.34) yields the following theorem, which is simply a rephrasing of an homolo- 
gous theorem in [7], with the proper extension of understanding to the setting of RHPs. 

Theorem 2.1 (Thm. 4.1 in [7]) Given two RHPs related by the elementary Schlesinger transformation 

(Def. 2.4), the difference of the Malgrange differential on the two solutions is a closed differential on the 

deformation manifold given by 

u M (d; [f]) - cu M (d; [T]) = dlnH jj J j (2.50) 
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where 



H 



i 3 



(T(»?))« 
(rtf)" 1 )* 



for (i^j) 
/or {7/}, r^oo 

/or e^co 
for {: ; ? } . x , 



(2.51) 



where the notation r'(oo) -as previously- denotes the derivative in the local parameter 

T(z) =: l + -r'(oo) + ... 



(2.52) 



For the reader's convenience we verify Thm. 2.1 for the case { $ J } with distinct £ ^ oo ^= r\ and both £, r\ £ £7. 
Let d be a variation of the jump-matrices M (i.e. not moving £,77). We have from Prop. 2.3 



R -i R , = mE^T-^o 



(r-i(f)r(»?))ii V«-»? 

We then have to compute the (possibly formal) residue (2.34) 



(2.53) 



res Tr 



9r(z)r- 1 (z) 



(r^gr^jg (r- 1 (par(e)r- 1 (g)r( ?? )) i 
(r-i(Or(ry)k 



(r-^orfa)), 



(r-i(Oar(»?))« , (5(r-H£))r(r?)) i 



(r-HOr(v?))i 



9 In 



(r-KOrW)ii 
(r-^or^jy 



£-?7 



(2.54) 



To verify the formula for d^, as well we must add a small circle around them and a new jump M(z) = (z — £) Ei , 
M(z) = (z — rj)~ Ej respectively. Since £,77 did not exist as deformation parameters in the RHP for V, from the 
definition of uj m we need to compute (we do it only for d^, leaving the verification for to the reader) 



r (%[r]) = - 



Tr I ri 1 ^ 



E, 



dx 



X — £ / 2z7T 



(2.55) 



l\z=t\=e 

Note that f _ = RT (and f + = RT (z - 0~ Ej ) so that 

'R- 1 {z)R'{z)T{z)E l T-\z) , T-\z)T'(z)E i 



w M (%[r]) = -resTr 



1 



(r- 1 (£)r'(£)r- 1 (£)r(r ? )) i 



d( In 



'(r-HOrfo)) 



This proves completely the case considered. It appears quite obviously that 



(2.56) 
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Proposition 2.4 An elementary Schlesinger transformation exists if and only ifH^j X =/= 0. 

Remark 2.8 The zeroes of the matrix entries of the solution of the RHP acquire therefore the meaning of in- 
tersection of the Malgrange divisor with space of the parameter added to the problem, namely, the position of the 
Schlesinger transform I °? V X . 

We are not going to dwell at length on the algebra of iterated elementary Schlesinger transformations and on the 
general transformation since the formulae are contained in [7]; we only point out that in loc. cit. the transformations 
were applied to solution of either isomonodromic or isospectral deformation problems and not a general Ricmann- 
Hilbcrt problem. Thus, we are mainly shifting the perspective (and compactifying some notation) of [7]. 

Remark 2.9 The fact that to evaluated on two solution of RHPs with the same jumps is a closed differential is 
immediate from the fact that the curvature of lo does not depend on the growth behavior of the solution T but 
only on the jump matrices. When we specialize the setting to the case relevant to ODEs, Thm. 2.1 will hold for 
differentiations with respect the monodromy data as well. 

2.3.1 Generalized Sato formula 

Let r\ £ £7: denote by T rj . i j(z) the Schlesinger transform 1°°^' j of the solution of the RHP 2.1. Then the second 
formula (2.51) reads 

r«(r?) ocexp / w M (»;[r„ ;i j]) - ui M (»;[T]) (2.57) 

where the one form under integration is closed by the above remark. This is nothing but Sato formula for the Baker 
Akhiezer vector; of course, at this level of generality we do not have a "tau" function because -in general- lo m will 
not be a closed differential of the deformation parameters of the problem. If the problem admits a r-function, 
that is, if the differential lo m (or f2) is closed on the submanifold of Ricmann Hilbcrt problems under consideration 
then we have a honest version of Sato formula 

J OO T} 1 

r ij (r 1 ) = -^± (2.58) 

where r j °? V j stands for the r-function of the problem with the "insertion" of the Schlesinger transform. 

To make the remark a bit more concrete, let us pick a point z = a ^ £7 and a small disk D. For simplicity 
in writing the formulae we will simply set a = 0. Let T(z) be a diagonal matrix of the defined on <9B that admits 
analytic continuation on P 1 \ D 

00 „ 

T(z) = J2 -r ' Tk = dia gfei' • • • . *fc,r) . (2-59) 
fe=i 

where the Laurent series is supposed to be actually convergent on P 1 \ D. Let T(z; [T]) denote the solution of the 
RHP 

P+ = r_M , z e £7 (2.60) 

r+ =r_e T , zedo (2.61) 

r(oo) = 1 (2.62) 
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If the L°° n L 2 norm of e T — 1 is sufficiently small, then the solvability of the problem is guaranteed by standard 
perturbation theorems. Therefore T(z; [T]) is defined at least in a ball around T = 0. Since the jump matrix e T is 
diagonal (and analytic in the complement of the disk), it is immediately seen from Prop. 2.2 that lo m — f2 5 and 
that they are closed as differentials on the (infinite dimensional) manifold of T's. Thus there is a locally defined 
function such that 

5WT) := I Tr (n^T) — (2.63) 

J dB 2l7T 

Denote now by T | i j} (^) 6 the tau -function resulting after the elementary Schlesinger transformation j ^ V\. 
Then the content of Thm. 2.1 can be rephrased as 



r(T) 



f(0 



{ c ,7}( r ) 



r(T) 



(2.64) 



In fact, more is true: if £ falls within the disk D then the RHP for the Schlesinger transforms | °° ^ j , | | can 
be formulated as (exercise) 



r+ = r_c J i - 



Ml 



,Ei-Ei 



z&dB 



r + = r_r*'«M/< z e s 7 
r(oo) = l 



r + = r_c'(i-i)"'= 



r + = r_« B iMr^ zes 7 
r(oo) = l 



The jumps on 9D can be written 

-Ei 



1 



c 



exp 



E 

k=l 



T k + E 3 C k /k 



exp 



E 

fe=i 



This leads to the following identities 



r«(C; [T]) = r y (0; [t - ^[C]]) , r~VC; PI) = r tj (o ; [t + E&]]) 



(2.65) 



(2.66) 



(2.67) 



and hence 



r«(C; PI) 



-{°°/}(r-^-[c]) 

r(T) 



>2 

Here we have used the standard notation [Q — (£, ^- 



r«(C; PI) 



(2.68) 



2.3.2 Hirota bilinear relations 

We will not go into much depth here, since all is well-known but it may give an analytic perspective on the relations, 
which are usually taken only formally. 

As in the previous Section 2.3.1 let C be a counterclockwise circle around a point z = a £ £7 (a = for 
simplicity) and let T(z) : C -> gl r (C) be as (2.59). We have 



r + (C; [T])e-no+no T -\ C;[ f]) 



c 



2i7rC 2 



r_(C;[T])r: 1 (C; PW 



dC 



c 



2iTT( 2 



(2.69) 



5 This follows by observing that i9 will be identically zero due to the analyticity of T outside of the disk. 
6 The notation T stands for the (infinite) vector of the (matrix) coefficients of T(z), T = (Ti, T2, . . .). 
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Now, the matrix T(x; [T])T- 1 (x; [T]) has no jumps outside C since the other jumps of the problem have been 
left unmodified; thus it is analytic on the complement of the disk D and goes to the identity at oo. Thus we have 
the identity 

/ r + ( C ; [T])e-nO + nO T -^ i^j dC = Q (2 J0) 

Jc 2l7r C 

which is valid identically in T, T in a neighborhood of T = 0. 

In view of the interpretation of T as the Baker-Akhiezer vector for the r-function (Sato-formula above), the 

reader may regard (2.70) as the generating function of an infinity of bilinear identities between matrix-valued t 

functions; 



E 



{ " / } (T - EtlQ&O-Wr {° e j}(f + EM)f - (2-71) 



The identity (2.71) should be used as a generating function of an infinite hierarchy of PDEs for the matrix-valued 
tau-function [T(T)]ij — r | c fj ) | (T) when expanding it in Taylor series with respect to T around the diagonal 

T = T. This generates a sort of "addition theorem" for tau-functions (see Remark 2 in [8]). 

The variational formula (2.9) in Lemma 2.1 takes on an added significance in view of the identity 

1 



— T 



{!;}= 


[r-Hflrfo)] 







(2.72) 

ij 



To explain this in a simple situation we now consider two disks Do and Dj centered at two points -say- a = 0, 1; 
on the boundaries of these disks we introduce diagonal jumps exactly as in Sec. 2.3.1 

OO rp(0) OO rp(l) 

fe=l k=l ^ ' 

We denote by t(T (0) ,TW) the tau-function as a function of the two (infinite) sets of times and t j ° j | (f^\f^) 
the Schlcsinger-transformed one. Let £ € Dq and 77 € ©i; retracing the steps that lead to (2.68) we find 



( T (o) jT (i)) 



T 



(2.74) 



This formula is the content (in different notation) of (3.11) in Thm. 3.2 of [8]. 

The variational equation (2.9) in our Lemma 2.1 applied to derivatives with respect to some directions in T^°\ 
then become the generating functions for the Hirota bilinear relations that appear in Thm. 3.4 of [8]. They all 
boil down to the following identity, to be understood as generating functions of PDEs when evaluating its Taylor 
expansion on the diagonal T^> = j = 0, 1 

T- 1 (^T)T^(x;T)TZ 1 (x;f)T( m f) dx = 
£ — x x — 77 2iir 

L- 1 (e;T)L+(a;;r) ? m _ T m ^(x; f )T( m f ) dx £ L _1 (£; T)T + (x; T) f(.i)_ T {i)Y^ L {x;T)Y{ l q;T) dx 



-e 



^2.75) 



£-x x - r] 2in J d0i £-x x-t] 2i 

The interested reader should compare this with Theorem 3.4 in [8]. Since it is not the primary focus of this paper 
(and it is certainly not a new result), we will not pursue the issue here, also because it has been dealt with at 
length in [8] , even though in the context of isomonodromic and isospectral deformations only. 
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2.4 Right gauge equivalence 

We will say that the two problems are (right) gauge equivalent if there exists an analytic function 

G : C\£ 7 -> GL r (C) (2.76) 

admitting boundary values (also for its derivative) at £7 and such that the jump matrices stand in the relation 

M{z) = GZ 1 {z)M{z)G + {z) , z G £7. (2.77) 

It is immediate then that the two solutions are related by T(z) = T(z)G(z). It is then seen that the difference of 
u) M (or fl) along the two solution T, T differ only in terms that do not involve r or T and depend only and explicitly 
on M, G. Thus this equivalence will not modify the singularity locus of ui M and -if both deformation families admit 
a tau function- both tau functions will differ only by multiplication by a smooth nonzero factor. Indeed a direct 
computation yields 

u, M (d;[T])-cj M (d;[T})= / Tr (r^r'.Gl^GL - T^T' + GZ l dG + ) — + (2.78) 

+ / Tr (M~ 1 M'GZ 1 dG + - GZ l G'_GZ l dG_ - GZ 1 G'_dMM- 1 + GZ^GZMG^dG+M- 1 ) (2.79) 

Js 7 2ztt 

The term in (2.78) vanishes by Cauchy theorem, since it amounts to the (boundary value of) the integral of 
Tr(r _1 r"G _1 9G) on a collection of contours surrounding £7 and contractible. Thus the difference will not involve 
the solution of the RHP and hence be a smooth differential in the parameters (i.e. gauge equivalence cannot modify 
the Malgrange divisor). Such equivalence does play a role in some cases (see [1] for examples where this happens, 
although not phrased in these terms). 

3 Rational differential equations in terms of Riemann— Hilbert data 

We now describe the class of ODE's with rational coefficients 

= A(z)V(z) (3.1) 

in an unconventional way: we will start from the formulation of a RHP and then indicate how this problem (when 
solvable) is equivalent to an ODE. The class of matrices A(z) that will be eventually described has poles at points 
ax, • • • , dk with orders n\ + 1, . . . , + 1. If rij = then the pole is simple. 

One may take the point of view that we are providing a different (transcendental) coordinate system on the 
finite -dimensional vector space of rational matrices with fixed polar divisor. 

The forward problem, namely, the construction of the RHP from the matrix A(z) is more standard and we 
only sketch the main points, since it does not really play a direct role here. This procedure is often called the 
(extended) monodromy map. The standard reference for many assertions below is Wasow's book [17] but also 
the paper [9] provides a concise recall. 
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3.1 Monodromy map 

Given a rational matrix 

K nj+1 no-1 

A ^ -EE Tf^V + E ***** ( 3 - 2 ) 

we consider the ODE $'(z) = A(z)^(z). Without loss of generality we will assume TrA(z) = so that any solution 
has constant determinant that we can assume to be unity deff = 1. We make the usual assumption that 

Assumption 3.1 (Genericity) The leading- coefficient matrices Aj^ nj have distinct simple eigenvalues. The di- 
agonal matrix of eigenvalues (with an arbitrarily chosen order) will be denoted by Tj^ nj . In addition, if rij = 
(simple pole) then Tj }0 has eigenvalues which are also distinct modulo 7L (i.e. no pair of eigenvalues differ by an 
integer). 

To simplify some issues in the general description we will assume that in fact there is no pole of the connection 
d z — A(z) at z = oo; this can always be achieved without loss of generality by a Mobius transformation that maps 
ootoa finite point (without mapping any of the other poles to infinity!). This allows to choose as basepoint z the 
point at infinity. 

Monodromy representation. We choose a basepoint z or the homotopy group and consider the initial value 
problem 'I'(oo) = 1. By analytic continuation of the solution around a loop that "goes around a/' (i.e. has index 
one relative to aj and zero relative to all other poles) we obtain ^(z) i— » ^{z)MJ l , det Mj = 1. These loops 
generate the fundamental group 7Ti(C \ {a\, . . .}, zq) and provide a representation of this fundamental group 

n(P l \{a 1 ,...},zo)~SL r (C). (3.3) 

Note that the basepoint for the normalization (oo in our case) and the basepoint for the homotopy group may 
not be the same. We will denoted by V the simply connected domain of P 1 obtained by dissecting P 1 along 
nonintersecting smooth arcs joining zo with each of the poles. 

Stokes' phenomenon in brief. Consider a higher order pole aj, with rij > 1 and denote as follows the distinct 
eigenvalues of the leading coefficient matrix Aj >nj 

T jtnj = diag(Ai, . . . , A r ) , A, ^ Aj, j (3.4) 
C := (z — aj) , n := rij , (3.5) 

(or ( = i for the pole at infinity). One can find at each pole aj 2rij directions (the anti-Stokes directions) 
angularly separated by ^ such that along each o f them there exists a permutation a (uniquely defined) yielding 
the definite ordering below 

Note that if on a direction we have the above ordering, on the next (counter) clockwise we have the exact reversed 
(with the same permutation). 
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The theorem which can be found in several place (e.g. [17]) is the following; given our IVP *5f(z) on the simply 
connected domain T>, there are 2n matrices S v , v = 1 . . . In (the Stokes' matrices) which are of the form 

S v = P„(l + N^P- 1 (3.7) 

where N v is upper-triangular on the directions where the ordering is as in (3.6) and lower-triangular on the directions 
where the ordering is reversed, and P a is the permutation matrix corresponding to the permutation a appearing 
in (3.6). There is an invertible matrix C = Cj € SL r (C) (the connection matrix) and diagonal matrix 

n 

H(Q := e T ^( L , T(C) = T ^~ e ' C = ( z - «) (° r C = \ for the pole at oo) (3.8) 

e=i 

in which T n = Tj n . is the matrix of eigenvalues of the leading coefficient used in the definition of the anti-Stokes 
directions. This matrix will be called the toral element (for lack of a better naming, since it belongs to the 
complex toral subalgebra of SL r (C)) and the matrix L is called the exponents of formal monodromy. 

These matrices (Stokes' connections and toral elements) are uniquely determined by the following set of condi- 
tions 

1. In one (arbitrarily chosen and then fixed) sector separated by the two consecutive anti-Stokes directions we 
have the asymptotic expansion 

~ f(z)e T(c) C L C (3.9) 
Y(z) = G^ + ^Y j ( J j,detG = l (3.10) 

2. In the next sector counterclockwise we have 

9(z) ~ Y{z)e T ^C L S 1 C (3.11) 
where we have labeled by 1 the anti-Stokes direction separating the two sectors, and in general 

tf (*) ~ ?(z)e T ^C L ^ ■■■S 2 -S 1 C (3.12) 
where we are in the sector between two anti-Stokes and we have crossed v such lines. 

3. The monodromy around the pole under consideration, the connection matrix and the Stokes matrices satisfy 
the condition 

M = C- 1 S 2n ---S 1 C (3.13) 

Note that at a simple pole we simply have n = and hence there are no anti-Stokes' lines but only the exponents 
of formal monodromy L and the connection matrix C. 

Remark 3.1 (Notational issue) It should be understood that all the above matrices (C, S,T, L) are different 
for each of the poles and hence we should understand an index distinguishing them and related to the pole under 
consideration. 

Remark 3.2 (Important) The triangularity condition for the Stokes' matrices S u can be expressed intrinsically 
by saying that 

M u {z) := e T(z) ^e" T(2) = 1 + 0((z - a)°°) (3.14) 
as z approaches a along the corresponding anti-Stokes direction. 

Since these are standard facts about ODEs we will not dwell on other subtleties. 
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3.1.1 Forward Birkhoff map. 

Given the rational matrix A(z) of our form, a chosen basepoint and dissection T>, anti-Stokes lines and base-sector 
from where to start the counting we have associated the collection of all the data 

M '.= {a j ,T j {Cj),L j ,C j ,{S vd v = 1...2n j },} j=1 ,..., k (3.15) 

The statement is that the map is (locally) injective but not surjective in general; there are some choices of data 
in A4 for which there is no corresponding ODE of the specified form. As the reader may suspect (or know), these 
data constitute the Malgrange divisor. 

There are two logically (and historically) distinct types of data in the above: the isomonodromic times 

i ■■ K-'/.'A,'h, (3-16) 

(the positions aj and the coefficients of the matrices Tj(£)). On the other hand there arc the genuine generalized 
monodromy data 

S := {Lj,Cj, {S v<j v = 1 . . . 2n,-}, } 3 =i... k (3.17) 

and M. is locally the product of the two. Of course this separation can be only done locally since the anti-Stokes 
directions depend on the leading coefficients of Tj (() and the dissection V depends on the position of the poles, and 
hence some care is in order -if global question are at issue- in describing the patching of these local descriptions 
(see for example [5]). 

4 Inverse Birkhoff map: Riemann— Hilbert problem 

Given the Birkhoff data described above, the question arises as to whether one can invert the map: starting from a 
concrete dissection, anti-Stokes lines etc, together with all the matrices appearing in M, can one reconstruct A{z)l . 
To this end it is necessary to specify a Riemann-Hilbert problem in the same spirit as Section 2. 

4.1 The set of contours 

The set of contours E7 consists of (see Fig. 3) 

1. for each pole we draw a circle not containing any other pole: we will call this the connection circle; 

2. for each pole aj a smaller circle is chosen, called the formal monodromy circle. On this circle a point (3j 
is chosen. 

3. each j3j on the formal monodromy circle connected with a set of mutually nonintersecting paths (stems) to 
the basepoint; 

4. at the higher poles we choose a third smaller circle called the toral circle; 

5. a point (arbitrarily chosen) on the toral circle is connected finally to the pole by In smooth curves that 
approach the singularity along the directions mentioned in Section 3.1; 

The word description is awkward but Fig. 3 should clarify all the elements. 
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Figure 3: The typical set of contours. The dashed lines within the formal monodromy circles are not jumps, but 
the domain where the determination of the powers (z — a)~ L is continuously defined so as to give a precise meaning 
to the expressions for the various jump matrices. 

4.2 The jump-matrices 

Rather than describing the matrices in words, we refer to the picture (Fig. 3) where we depict a situation with 
only two points ax, a,2 with n\ = (simple pole). The general picture is quite simply a repetition of several copies 
of the basic elements already manifest here. Since there is no monodromy around the basepoint zq = oo we have 
to impose the constraint 

c -i e 2M n(Jn . . . c -i e 2M lCi = 1 ( 4J ) 

It will be understood that the matrices (z — a) L are defined on the formal monodromy circles as continuous functions 
taken away the point of insertion of the points (3j (Fig. 3). 

4.3 Riemann— Hilbert problem 

Problem 4.1 Find a piecewise analytic matrix-valued function T(z) on the complement of the indicated contours 
as per Fig. 3 so that 

• on each arc it solves T + (z) = r_(z)A/(z) with the matrix M{z) as indicated in Fig. 3; 

• it is bounded on P ; 

• satisfies the normalization r(oo) = 1. 
Some properties follow immediately: 
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• Any solution satisfies detr(z) = 1. Indeed detM(z) = 1 implies that det(r(^)) has no jumps across the 
contours. Since T(z) is bounded so must be detr and hence it is an entire function, bounded everywhere, 
hence a constant. Since T(zq) = 1 then detT(z) = detr(zo) = 1. 

• If a solution exists, it is unique: indeed any two solutions T, T are analytically invertible by the above remark. 
Hence IT -1 is promptly seen to have no jumps, be uniformly bounded, hence a constant. The normalization 
forces r = r. 

The connection with the ODE is as follows; define the piecewise 
analytic matrix ^(z) related to the solution V as follows: 

• ^f(z) — r(z) outside the formal monodromy circles; 

• ^(z) — T(z)(z — CLj) Lj in the annulus between toral and 
formal- monodromy circles; 

• ^(z) = T{z)e T i( zS> (z — aj) Lj inside the formal-monodromy 
circles. 

The cut of the logarithm is taken where the dotted line is traced 
in Fig. 3. It is promptly seen that ^(z) solves a new RHP where 
the jump matrices are piecewise constant and unimodular. Thus 




Figure 4: The constant-jump Riemann- 
Hilbert problem for "P. 



A(z) := ^'(z)*- 1 ^) 



(4.2) 



is a meromorphic (traceless) matrix function with isolated singularities at the Oj's. A local inspection shows that 
it has poles of finite order there and hence it is rational. 



5 Tau function and (iso)monodromic deformations 

The goal now is to show that the modified Malgrange form is closed when evaluated along the manifold C C Q 
consisting of jump-matrices M(z)'s described in the previous section. We emphasize that the parameters are 

• the toral data Tj(z)'s; 

• the connection matrices Cj\ 

• the Stokes' matrices S v f : 

• the exponents of formal monodromy Lj . 

The reader acquainted with the literature about isomonodromic deformations should now realize that we allow 
many more directions of deformations. The computation of the closure of f2 relies directly on Prop. 2.2: 

SQ(d, 8) = dn(d) - tkl(d) = \ I Tr (m'm- 1 

It is clear that no contribution to (5.1) can come from the arcs where M(z) is independent of z. This leaves only 
the contributions coming from the toral circles, the formal monodromy circles and the Stokes' lines. Let us consider 
each type separately 



dMM~\dMM- 



ax 

2iTT 



(5.1) 
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• Stokes' lines; up to a constant conjugation by a permutation matrix the jumps are of the form 1 + N(z) 
with a strictly triangular N(z), and hence they do not contribute to (5.1) (the trace is identically zero); 

• Toral circles; on the Toral circles the jumps are diagonal and hence the commutator in (5.1) is identically 
zero. 

• Formal monodromy circles; once more, since the jumps are diagonal, no contribution to (5.1) comes from 
them. Note that at a simple pole the matrix L may actually be upper triangular: in this case the trace 
vanishes identically. 

Therefore we have the 

Theorem 5.1 The differential restricted to the submanifold C of group -valued jump -matrices described above is 
closed and defines a local function via the formula 

T(T,a,L,S,C) = cf n (5.2) 

This function is defined up to nonzero multiplicative constant and it vanishes precisely when the Riemann-Hilbert 
problem is not solvable, namely, on the Malgrange divisor. 

The fact that r has zeroes (and not a branching singularity) does not follow from our construction, it follows from 
[14] and [11]. 

We would like to stress that, had we chosen the Malgrange form u M directly, the closure would have failed. 
Indeed the reader may check that 

Sco M (d, 0) = ~ £ (da 3 dTv{L*) - da 3 dTv{L])) K 3 (5.3) 
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hx{x — a,j) dx 1 



where the integrals are on the toral circles with basepoint at the necks j3j. 

Aside from this, the correction therm •& in (2.25) does not contribute anything for most of the contours in £7; in 
fact the only contribution comes -not surprisingly in view of the above computation- from the formal monodromy 
circles 



1 f , 1 „ 1 x dx 1 v-^ , „ , f ln(x — a) dx , 
i?(9) = -/ TriM'M^dMM- 1 )— = - V Tr (LdL ) <b — '-— = (5.5) 

2 J-£ 7 2iT: 2 ^—^ J p. x — a 2in 

= \Y. Tr(LdL)(2in + 21n(& - a,)) (5.6) 
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Note that this "problem" is invisible if we allow only isomonodromic deformations, in which case using oj m or 
yields the same differential. 

Before we turn to a list of applications of the above result to explicit examples, we make the connection with 
the definition of [9] . 

5.1 Relationship between uj and uj 
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In order to establish the relationship we must "freeze" all the mon- 
odromy part of the Birkhoff data, that is, the connection matrices, 
the Stokes matrices and the exponents of formal monodromy. As 
noted above, in this case f2 = lo m . Then the observation that u M 
and uj JMU coincide can be dug out [14] but we re-derive it here 
for the sake of self-containcdncss. If now d is a derivative along 
the isomonodromic submanifold of the Birkhoff data, the differ- 
ential reduces to an integral only on the Stokes' lines, the toral 
circles and the formal- monodromy circles. Since the expression is 
repeated for each pole, we consider only the localization at one pole 
a = a,j, i.e. we vary the toral data/position only of one pole. Let d 
denote one such deformation involving only data at a — a,j . Then 

da; 



fi(0) = uj M {d) = - J Tt^Z l V_dMM- 



2i?r 



(5.7) 




There are two types of integrals here: the integrals along the 
Stokes lines and the integral around the toral circle. Along each 
Stokes contour the integrand can be equivalently written as follows 



Figure 5: The set of contours for the compu- 
tation of U) ..... 



Tr 



r -i r / 



M 



df - e T 5e- T df & T S~ 1 q- t 



J J 



V V 

= Tr ( I I ' df - I . 1 1 '. i), j \ i^TdT) 

T(z) := T(z) + Lln(z- a) 
where we have used the cyclicity of the trace on the second term followed by 

M _1 r~ 1 r / m = rv 1 r / , - m^m' 



0(t 



(5.8) 
(5.9) 

(5.10) 



Recall also that expressions like dT do not contain a logarithmic term because dL = (only isomonodromic 
deformations are allowed). As a consequence of (5.8) the resulting sum of integrals (refer to Fig. 5) can be 
evaluated by excising an e circle 7 around the pole z — a and integrating Tr(r -1 r'<9T) on the solid black contours 
indicated in Fig. 5, followed by a limit e —> 0. For fixed e the integral along the solid "tentacles" is equal (by 
Cauchy's theorem) to the counterclockwise integral along the e circle and the clockwise integral along the bigger 
dashed circle indicated in Fig. 5. 

The two integrals of Tr(r _1 r'9T) along the two circles surrounding the toral circle reduce -due to the jump of 
r there to 

,0 



Tr(rz TLdT - TT L T' + dT) dx 



toral 



dT) 



'toral 



& Tr 

J toral 



r -i r , 



Lda 



dx 
2in 



(5.11) 



7 Note that T has a jump of the form e^^ L in the straight ray from a to the toral circle (Fig. 3), but the expression Tr(r 1 r'<9T) 
does not have a jump there, since L and T are both diagonal. This is the reason why we did not draw the corresponding solid line 
in Fig. 5. 
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where the first integral is zero because it is residueless. The remaining can be deformed within the annulus between 
the toral and formal-monodromy circle (and the T_ of the toral circle becomes the r + of the local monodromy 
circle) and combines with the outer integral therein to give 

Lda \ dx f Tr(L 2 )da 2 dx 

(z — a) 2 2iir 



for.mon. 



Tr T^T' 



r -i r / 



a — x I 2iir 







(5.12) 



for.mon. 



The remaining integrals can therefore be retracted to an integral around the e circle: 

da; 



Tr T~ 1 T'dT 



2iir 



(5.13) 



' \z— a\=e 

In each sector we have the expansion valid for any N 

T(z) = G (l + £ - a) k ^ + 0((z - a) N+1 ) ^ll + f] Y k (z - a) k ^ (5.14) 

where the coefficient matrices G, are the same irrespectively of the sector. If we replace T(z) in (5.13) by a 
suitably high truncation of the formal scries we commit an error that it is easily estimated to vanish as e tends to 
zero. On the other hand the new integral is independent of e and reduces to the formal residue 



/ oo 

n(d; [T]) = u M (d- [V]) = res Tr (y" 1 (z)Y' (z)dT (z)) , Y := G 1 + £y fc ( z _ a) 1 



(5.15) 



to be understood simply as the coefficient of (z — a) -1 of the above formal series with a finite Laurent tail. The 
reader acquainted with the definition of lojmu will recognize that the expression is precisely the same defining 
u) JUU in [9] (of course one should repeat the residue computation at each pole and sum up). 
In conclusion we have shown 

Proposition 5.1 The (modified) Malgrange differential f2 restricted to the manifold of isomonodromic deforma- 
tions coincides with the Jimbo-Miwa-Ueno differential. 



6 Applications and examples 
6.1 Painleve II equation 

We single out the second Painleve for its relative simplicity in formulas and as an illustration of potential applications 
to other Painleve equations. We follow the setup in [6]. The Riemann-Hilbert description is depicted in the picture, 
with the notations and condition 



loo 



1 







se fz 3 +itz x 



U(s) 



1 se a 



1 j ^ 

s% - s 2 + s 3 + sis 2 s 3 = (6.2) 

Condition 6.2 is the condition that the products of the jumps at the origin is the identity. The Riemann-Hilbert 
problem is then that of finding T(z) uniformly bounded, with the indicated jumps and the normalization condition 

r(z) - l + oiz- 1 ) 

The matrix T(z) and the Stokes' matrices are simply 



(6.3) 



T ( z ) = ~ [ tt 2 - 3 + itz ) 03 , 03 := 



1 

-1 



(6.4) 
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and the matrix ^(z) := T(z)e T has constant jumps (on the 
anti-Stokes lines) and satisfies 



*'(z)*(z) _1 = -i (iz 2 + t + 2u 2 ) cr 3 - 4uzer 2 - 2va x 



-U 

o a 



z 2 + 



4m 
—Aiu 



z + 





' 


1 " 




' -i~ 


a x := 


1 





, cr 2 := 


i 



-it — 2iu 2 —2v , 

-2v tt + Ku 2 H 6 - 5 ^ 



, « = </(*)(6.6) 



Under an isomonodromic deformation the coefficients u, v in (6.5) 
must evolve as functions of x and -in particular- 



u = u(t; s) = 2 lim zT\2(z; i, s) 



(6.7) 



(the limit does not depend on the sector we choose) solves the 
second Painleve equation 



d 2 u 
dt 2 



2u J + tu 



(6.8) 




By direct computations we have 



duY 
~dl ) 



tu 2 - u 4 dt = d t In r(t; s) dt , d 2 In r(i; s) = u(t; sf 



(6.9) 



Gcnerically we can solve the condition (6.2) for S3 and use s\, s 2 as independent variables: more appropriately one 
should consider the parameter space as the algebraic manifold specified by (6.2). Introduce the matrix kernel 

r(x)- 1 r( y ) 



K(x,y) ■, 



x-y 



(6.10) 



We then have 
d Sl lnr(i; s) 



Tr r _1 r 



i4 3 1 -j. 

-^-X -j-ltX Q 



1 + sl 



dx 



2iir (1 + S1S2) 2 Ji 



Tr r _1 r 



q — itx 



dx 
2in 



= - Ki 2 (x,a;)e 



1 + si 



dx 



2in (1 + sis 2 ) 2 Ji 



K 2 i(a;, x)e~ 



, dx 

2iTT 



(6.11) 



where the boundary-value indication is irrelevant since the indicated matrix element does not have a jump on 
the corresponding line. Using formula (2.15) for the second derivatives and then (2.9) one may derive an integral 
formula for the derivative <9 2 9 Sl lnT; however, in view of (6.7) it is simpler to use the variational formula (2.8) 
directly, thus yielding 



2 lim 



z r(x) 



d Sl u(t;s) = 2 lim zd Sl Ti 2 

z — >oo 





■^x 3 -\-itx Q 



T-\x)T{z) 



12 



da; 



-2 lim 



(l + sis 2 ) 



z T(x) 



x — z 
e — ^ x 




2in 



T-^x)T(z) 



12 



da; 



x — z 



2in 



(6.12) 
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Since T(z) — » 1 we have 

d Sl u(t;s) 



(I + S1S2) 2 J 1 



r n (a;) 2 e-^ 



1. -^^-itx d% 2 



2i7T 



r 12 ( a; ) 2 eT- 3 +^ — 



(6.13) 



6.2 Variation of finite Toeplitz determinants for discontinuous symbols 

Let S 1 — {\z\ = 1}. Suppose that we have the following data: 

• A finite partition of S 1 in subarcs 7j, with (ij being the separating points; 

• a collection of functions a,j(z) : 7^ — > C which are analytic in a neighborhood of "fj. 

Define a(z) : S 1 —> C as the piecewise analytic function that coincides with ctj on each (interior of) jj. The n-th 
Toeplitz determinant is defined as 

dz 



T n [a] := det\p,j-i] it j =1 . 



t'j 



z J a(z) 



1*1=1 



2inz 



, T [o]:=l 



(6.14) 



It was shown in [1] that T n [a] (and more general objects) are isomonodromic tau functions for the following RHP; 



r» = r n (z) , r+ = r_ 



1 w(z) 
1 



T(z)=T n (z)= (l + G(z- 1 ))z" CT3 
T(z) = C(ln(z - /3j)) , to(«) := z^a^z) 



(6.15) 



In order to reduce to the general setting of Sect. 2 we should add a big counterclockwise circle (the "formal 
monodromy" circle) with jump z n<73 and replace the asymptotics with r(oo) = 1; to take care of the logarithmic 
growth at the endpoints f3j we take e-circles centered at the f3j (with e sufficiently small) and re-define V as follows 



T(z) :=T(z) 



1 -C[w](z) 
1 



, \z-fy\<e, C[w](z) := 



w(x) dx 
| = 1 x — z 2i?r 



(6.16) 



Thus we recast the problem (6.15) in the equivalent one 



r+ 


= r_ 


r+ 


= r_ 


r+ 





1 z n a(z) 

1 

1 -C[w](z) 
1 

, \z\ = R > 1 

T(z)=T n (z)= (l + Oiz- 1 )) 



, \z\ = 1, min \z — j3j\ > e 
3 

\x-0j\=e 



(6.17) 




and r(z) uniformly bounded in C. 

Note that the additional jump on \z\ = R (oriented counter- 
clockwise) is independent of the symbol and hence undergoes no 
deformation. Therefore this has no impact in the definition of the 
differential fi„(«) = [T n ]). 

Note that the first column of the solution T to problem (6.17) guarantee that the solution is bounded every- 
consists of polynomials (of degree < n) for Izl < R; indeed it wnere - here C[w](2;) .— § x _ z 2ijr and the new 

RHP has no jump within the small disk cen- 
tered at [3j. 
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Figure 6: The modification to the RHP to 



has no jumps on \z\ — 1 and the added small circles (due to the 
triangularity of the jump-matrices). Outside the circle \z\ = R. 
said column divided by z n is bounded at infinity, which proves the assertion. 

Similarly, the second row of (which solves a RHP with the jump matrix on the left) also consists of 
polynomials, by parallel arguments. 

Thus the integrand in the definition of the differential f2 is in fact the integral of a polynomial kernel 8 

K 2 i (z, w) := ± ' - li21 (6.18) 

z — w 

of degree < n — 1 in both variables (it has no jumps, it is regular on the diagonal z — w and bounded by a power 
at infinity separately in each variable). 

Now, the computations in [1] where done for an isomonodromic tau function, however the proof passes through 
without change to the extended ^-differential. Namely, the proof consisted in showing that 

fi n+1 -n„=<nn(^±iM) (6.19) 
V T n[a] J 

where fi„ denotes the differential evaluated on the solution T n of (6.17). Since the two RHPs differ by an elementary 
Schlesinger transformation of the type discussed in Sect. 2.3, this part is identical. Moreover, f^o = because the 
integrands in (2.25) and (2.7) are identically zero (all the matrices are upper triangular and Sg is strictly upper 
triangular), which means that To = exp J £Iq is a constant independent of the symbol, and can be taken 1. Since 
T [a] = 1 as well, this implies that 

8hiT n [a}=n n (6.20) 
where 6 is the total differential w.r.t. any parameters may appear in the symbol and the positions of the endpoints. 



Dependence on the f3j: the dependence on each of the 0j is only in the jumps around the corresponding circle, 
and we have 







-Aj(/3) 1 





(6.21) 



where the subscripts L R indicate the boundary value on the Left/Right of w along the contour, relative to the 
orientation of the contour, namely, A is the jump of the symbol. This shows -a well-known fact- that 



dp, T n [a] = - f K 2 i(x,x) 

J\x-Pj\=e 



-AjjP) 1 dx = AjQg) 
2i7r x — Pj 2in 2in 



Kzi^-,^) 



(6.22) 



Dependence on the symbol. Let d denote a deformation of one of the 
functions af, then the jump matrices M(z) deforms only on the arc "fj and the 

8 A fact which is well known and key to matrix model computations. 
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Figure 7: The "dumbbell" contour 



small circles we have added, on a "dumb-bell" contour (Fig. 7). In particular, 
on the two circles the integral for reduces to 

/ K 21 (x,x)C[dw](x)^- 

x n daj(x; s) dx 



C[dw](z) 



x — z 



2in 



(6.23) 
(6.24) 



Since, ~K-2i(x,x) is a analytic in the interior of the e-circle (in fact, it is a 
polynomial!) then we can perform the integration by "contour deformation" 

leaving an integral of J£.2i{x, x) against the jump of C[3w]. In conclusion, taking into consideration all the contours 
we have simply 

r r\x 
dT n [a] = / K 2 i{x,x)x n daAx;s)— (6.25) 

In general we have simply 

dx 



d\nT n [a] 



1*1=1 



Tr r-^aor'te) 



x n d s a(x; s) 




2in 



= <t K 2 i(x,x)x n d s a(x;s) — 



K(x,y) 



dx 

2iTT 

r-H^rfa) 

x-y 



(6.26) 
(6.27) 
(6.28) 



Example 6.1 A particular case is ifa,j(x;s) = Sj are constants, in which case the above reduces to an integral on the 
sub- arc jj 

d Sj \nT n [a] = f K 21 (x,x)x n ^ (6.29) 
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6.3 Hankel and shifted Toplitz determinants 

The notion of semiclassical symbols was defined in [12] and it was shown in [1] that the corresponding Han- 
kel/Toeplitz determinant were isomonodromic tau function, to within an explicit non-vanishing factor. 

The weights we are considering are all of the semiclassical type as defined in [10, 12, 13, 2, 4]. This means that 
they are of the form n(x) = e~ v ^ with V'(x) an arbitrary rational function. 

They are integrated over contours 7j which can be arbitrary contours in the complex plane as long as all 
integrals x k n(x) dx are convergent integrals. The range for k will be either N or Z, depending on the situation; 
a detailed description of the contours can be found in [2, 3] and we refer thereto for a more detailed discussion. 

We will choose arbitrary complex constants xj for each contour "fj and use the notation 



f*V(*)d*:=£*J*V(^ = ^ 



(6.30) 



We will also use the notation k : C — > C to indicate the locally constant function that takes the (constant) value 
Kj on the corresponding contour jj. The Markov function (sometimes referred to as Weyl function) for these 
semiclassical weights is simply defined as the locally analytic function on C \ Ujj given by 

M(C)dC 



W(x) 



(6.31) 
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The function W(x) has logarithmic growth at the hard— edges i.e. the endpoints of contours jj where fi is 0(1). 
In this case one verifies that W(x) = C(ln \x — et|), where a is the hard-edge point. 
Consider the jump matrix supported on the contours (J 7^ 



M{z) 



1 



z e 7r 



1 

Consider the following Riemann-Hilbert problems for a matrix T = r*™'^ parametrized by two integers n, I 

r+ =r_M 



(6.32) 



r(z) = 



O(i) 

(l + ^z- 1 )) 
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z™-^ 1 

z-^ 1 



z->0 







(6.33) 



and r(z) = 0(hx{z — a)) at any "hard-edge". 

By definition of Schlesinger transformations (Sect. 2.3) we see that all these problems are successive Schlesinger 
transformations of the problem n — 0, t — 1, which has the immediate solution 



Y f c~ v,(a,) da; 
J x x—z 2in 

1 



(6.34) 



It is also apparent that the differential lo m (or f2, they coincide here because the term {) in Def. 2.2 vanishes 
identically, being the jumps 1+upper triangular) are closed differentials in the deformations of V, the endpoints 



and in the x/s. 



Moreover, comparing the change u> M (d; [T^ n ^]) between two problems with n^ii±lorl^f±l one sees 
exactly as in [1] that 

d\nAi=w M (d;[T^]) (6.35) 



where are the shifter Toeplitz determinants 



det 



/ H 



1 , A £ 







Here the case of Hankel determinants corresponds to I = n — 1; then A™ -1 is -up to a permutation of columns, 
hence a sign- a Hankel determinant. 

In particular the derivatives w.r.t. the parameters Xj are 



^.lnAl= f K^(x,x)e- V ^^- 



K^(x,y) = 



[rM(ic)-irM(»)] B 



x-y 



This kernel is the usual Christoffel-Darboux kernel: for example, if t = n — 1 (Hankel determinants) then 



a; - j/ 



j=0 J 



(6.36) 
(6.37) 

(6.38) 
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where pj are the monic orthogonal polynomials of exact degree j for the moment-functional, namely, 

Pj {x) Pk {x)e- v{x) dx = hjSjk (6.39) 

Remark 6.1 Note that the RHP for r*™^ can be converted to a RHP with constant jumps for ^ := re~^ CT3 and 
hence reduced to a rational ODE v]/'v|/ _1 = D(z). The resulting isomonodromic tau function is not exactly equal 
to A^ because the RHP for 1^™'^ and the RHP in the canonical form described in Sec. 4 differ by a gauge in the 
sense of Sect. 2.4- However the corresponding factor is easy to compute and this was accomplished in [1J. 

We conclude with the remark that the variational formula: for A^ are valid for very general weights, not necessarily 
of semiclassical type. Indeed we are really using only the results of Sect. 2, where the connection to rational ODE 
is not necessary. 
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